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A Spitzer Space Telescope Study of SN 2003gd: Still No Direct 
Evidence that Core-Collapse Supernovae are Major Dust Factories 
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ABSTRACT 

We present a new, detailed analysis of late-time mid-infrared (IR) observa- 
tions of the Type II-P supernova (SN) 2003gd. At about 16 months after the 
explosion, the mid-IR flux is consistent with emission from 4 x 10~^ M0 of newly 
condensed dust in the ejecta. At 22 months emission from point-like sources 
close to the SN position was detected at 8 //m and 24 /im. By 42 months the 
24 fim flux had faded. Considerations of luminosity and source size rule out the 
ejecta of SN 2003gd as the main origin of the emission at 22 months. A possible 
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alternative explanation for the emission at this later epoch is an IR echo from 
pre-exis ting circumstellar or inte rstellar dust. We conclude that, contrary to the 
claim of ISugerman et al. I (120061 ) . the mid-IR emission from SN 2003gd does not 
support the presence of 0.02 Mq of newly formed dust in the ejecta. There is, as 
yet, no direct evidence that core-collapse supernovae are major dust factories. 



Subject headings: supernovae: general — supernovae: individual (SN 2003gd) 



1. Introduction 



Massive stars explode via core collapse and ejection of their surrounding layers (e.g., 
Arnett et al. 1989, and references therein). The extent to which core-collapse super- 
novae (CCSNe) are, or have been, a major source of dust in the universe is of great in- 
terest. For many years it has been hypothesized that the physical conditions in the ex- 
panding ejecta of CCSNe could result in the condensation of large masses of dust grains 



pandmg e.iecta oi UUblNe could result m the condensation ol large masses ol dust 
( Cernuschi. Marsicano. fc Codina Ill967l: Hovle fc Wickramasinghe Ill9"70 : Gehrz 1^989 



199' 



Dwek 



19981: 



Todini fc Ferrara 



2001 



Nozawa et al. 



2003f ) 



Tielens 



This follows from the fact 
that large abundances of suitable refractory elements are present. In addition, cooling by 
adiabatic expansion and molecular emission takes place, and dynamical instabilities can pro- 
duce density enhancements or "clumping." This, in turn, will aid dust formation through 
the effects of cooling and self-shielding. Further support for these ideas is provided by iso- 
topic anomalies in me teorites which indicate that some grains must have formed in CCSNe 
jciavton et al. \\l99^ . 



Interest in CCSNe as dust producers has in creased recently due to the problem of ac- 

couiiting for the presence of dust a t high redshifts ( FalL Pei. &: McMahon Ill989 ; Fall. Chariot, fc Pei 



19961 : IPei. Fall, fc Bechtold Ill99ll : IPettini et al. Ill997 



Bertoldi et al. II2003I ). In these early 



eras, much less dust production from novae and asymptotic giant branch stars is expected 
since fewer stars will have evolved past the main-sequence phase. Consequently, CCSNe 
arisin g from Population III s t ars have been propos ed as the main early-universe source of 



dust (jTodini fc Ferrara 
(iTodini fc Ferrara,2001l : 



2001 



Nozawa et al. 



Nozawa et al. I l2003l ). Models of dust formation in CCSNe 



20031 ) succeed in producing copious amounts of dust 



around 0.1-1 Mq even in the low-metallicity environments at high redshifts. This corre- 



spon ds to a supernova (SN) dust condensation efficiency of about 0.2 (IMorgan fc Edmunds 



20031 ). where the efficiency is defined as the dust mass divided by the total mass of refrac- 



tory elements. This is enough to account for the quantity of dust seen at high redshifts (see 
Appendix). 
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Newly condensed dust in CCSNe can be detected by its attenuating effects on optical/near- 
infrared (IR) light or via thermal emission from the grains in the ejecta. These methods have 
been used in attempts to measure the dust productivity of CCSNe. Both methods are sub- 
ject to uncertainties due to dust formation in optically thick clumps, so the derived masses 
tend to be just lower limits. By far the most extensive evidence for ejecta dust condensation 



is that obtained from 


;he peculiar Type II Sr 


■vf 1987A, where both techniques were em- 


Dloved (Danzieer et al. 


19891: Lucv et al. 


1989: 


Meikle et al. 


19891: Iwhitelock et al. 19891: 


Suntzeff &; Bouchet 


1990l: Dwek et al. 


199 


2I: Roche. Aitken. & Smith 1993: 


Wooden et al. 


1993: Ercolano. Barlow. & Sueerman 


200 


7). However, even the highest value obtained is 


only 7.5 x 10"^ 


(Ercolano. Barlow. &; Sueerman 2007). 


Pozzo et al. ( 


2004) used the 



attenuation method to infer a dust mass exceeding 2 x 10"^ Mq in the Type Iln SN 1998S. 
However, such events are relatively rare. Moreover, in this case it is suggested that the dust 
condensation was not in the body of the ejecta, but rather took place in the cool, dense shell 
produced by the impact of the SN ejecta with circumstellar material (CSM). We note also 



that an alternatiy e IR-echo scenario for SN 1998S is not ruled out (iGerardy et al. 1 12002 



Pozzo et al. Il2004f ). 



Prior to the launch of the Spitzer Space Telescope (hereaft er, Spitzer), the o n ly ev i- 



dence of dust condensation in a typical CCSN was presented by lElmhamdi et al. I (120031 ) . 
who used optical line attenuation to infer a dust mass lower limit of about 10~^ Mq in the 
Type Il -plateau (H-P) SN 1999em. Mid-IR studies of the Cassiopeia A supernova remnant 
(SNR) JPwek et al. Ill987l : Eaeage et al. 1^9961 : IPouvion. Laeage. fc Pantin IboOlh indicate 
that dust formation took place during its explosion, but again the mass o f directly observed 



dust is small. Sub-millimeter studies of this SNR by iDunne et al. I ((20031) usi ng SCUBA led 



them to claim that at least 2 Mq of dust formed in the supernova. However, iKrause et al. 



(I2OO4J ) have used the same data together with observations from Spitzer to show that most 
of this emission originates fr om a line-of-sight molecular cloud, and not from dust formed in 
Cas A. iTemim et al. I (j2006[ ) used Spitzer observations to estimate 10~^ — 10~^ Mq of dust in 
the Crab Nebula SNR. While rather uncertain, this result may be more relevant to this paper 
than that of Ca s A, since the Crab Nebula is thought to have arisen from a progenitor o f 
mass 8-10 Mq (INomoto. Sugimoto. fc Sparks Ill982l : iKitaura. Janka. fc Hillebrandt 1 120061 ) . 
similar to that of the CCSN studied here (SN 2003gd). 



In summary, prior to the launch of Spitzer, direct observations of CCSNe or SNRs have 
never revealed more than ~10~^ Mq of dust — only ~1% of the mass required if CCSNe are 
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to be important dust sources. But the number of CCSNe investigated for dust production 
is small, and with the exception of SN 1999em, rather atypical. The availability of Spitzer 
has provided an excellent opportunity for us to test the ubiquity of dust condensation in a 
statistically significant number of typical CCSNe. It provides high-sensitivity imaging over 
the mid-IR, covering the likely peak of the dust thermal emission spectrum. This can provide 
a superior measure of the total flux, temperature and, possibly, dust emissivity than can be 
achieved at shorter wavelengths. Moreover, the longer-wavelength coverage of Spitzer lets 
us detect cooler grains and see more deeply into dust clumps than was previously possible 
for typical nearby CCSNe. In addition, multi-epoch observations with Spitzer can distin- 
guish between dust condensation and IR echoes via the strength and shape of the light curve. 



In this p aper we analyz e o bservations of the Type II- P SN 2003gd at three late- 

time epochs. iHendry et al. I (120051 ) and lSugerman et al. I (120061 ) (henceforth "S06") reported 
optical attenuation effects in the late-time spectra and BR light curves of SN 2003gd which 
indicate dust condensation in this event. Using Sj>zfeer observations at two late-time epochs, 
S06 also report mid-IR emission from the condensing dust. This was the first-ever report 
of condensing dust in a SN II-P on the basis of thermal emission from the grains. Here 
we present a new study of these Spitzer observations. We agree with S06 that some of the 
earlier-epoch mid-IR emission was due to a modest quantity of ejecta dust. However, we 
find that their principal conclusion, that the later-epoch observations indicate the presence 
of 0.02 Mq of dust formed in the ejecta, is not supported by the data. Consequently, thus 
far there is no direct evidence that CCSNe are major dust factories. 



Observations 



SN 2003gd was discovered (lEvans fc McNaught 1120031 ) on 2003 June 12 (UT dates are 
used throughout this paper) in t he SA(s)c galaxy NGC 6 28 (M74). On 2003 June 13 it 
was identified as a Type II event (iGarnavich &: Bass I l2003l ) using a J-band spectrum. On 
2003 June 14 the identification was confirmed using opt ical spectra, a nd it was estimated 
that the SN was roughly 1 month (jPhillips et al. 1120031 ) or 2 months (IKotak et al. II2003I ) 
post-explosion at t he time of discovery. Using light - curve comparison w ith other SNe II- 
P, it was deduced JVan Dvk. Li. fc Filippen"koll2003l : iHendrv et alTlboOsh that SN 2003gd 
was a normal Type II-P event with estirn ated explosi on dates of, respecti vely, 2003 March 



17 ± 3 (I Van Dvk. Li. fc Fihppenko Il2003f ) or 18 ± 21 (IHendrv et al. Il2005f ). We adopt 2003 



March 17 as the explosion date, 87 days pre- discovery. 
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On the basis of a variety of methods (standardized candle method, b rightest super- 



giant s, kinematic) a distance to SN 2003gd of 9.3 ± 1.8 Mpc was found (Hendry et al 



20051). ModeUing of the hght echo of SN 2003gd fIVan Dvk. Li. fc Fihppenko 



20061 ) suggests 



a somewhat smaller distance of about 7 Mpc. S06 adopted 9.3 Mpc, and so for ease of 
comparison with their work we shall adopt 9.3 Mpc thro ughout. Total extinction (Galactic 
+ host galaxy) estimates of E{B - V) = 0.13 ±0.03 mag JVan Dvk. Li. fc Filippeiiko]l2003h 
and E{B — V) = 0.14 ± 0.06 mag (IHendry et al. II2005I ) were reported. Using two indepen- 
dent methods (bolometric luminosity of exponential tail; direct comparison with SN 1987A 
bolometric light curve), Hendry et al. estimate an ej ected ^^Ni mass of 0.016ln'mK M^;^, only 
about a fifth of the ^^Ni mass found in SN 1987A Jwhitelock et al. 



1988; Bouchet et al. 



199ll ). The progenitor star was identified in archival images from the Hubble Space Telescope, 



the 2.6-m Nordic Optical Telescope, and Gemini North, as a red supergiant of mass 6-12 Mq 



fIVan Dvk. Li. fc Filippenko II2003I : ISmartt et al. 1120041 ) 



The field of SN 2003gd was observed with Spitzer's Infrared Array Camera (IRAC) at 
3.6, 4.5, 5.8, and 8.0 fim on 2004 July 25 and 28 (days 496 and 499) and again on 2005 
January 15 (day 670). The first two SN observations were obtaine d serendipitously withi n 



the Spitzer Infrared Nearby Galaxies Survey (SINGS) (PID: 0159) (IKennicutt et al. II2003I ). 



In the SINGS program observations are duphcated with a delay of a few d ays, to permit 



ident ification of asteroids and to better sample the emission on subpixel scales (IRegan et al. 



20041 ). In each wavelength channel the two images are combined to yield an "Enhanced 



Data Product," and these are publicly available from the NED database. The 2005 January 
observation was obtained within our Spitzer supernova program (PID: 3248). The SINGS 
program also used the Multiband Imaging Spectrometer for Spitzer (MIPS) to acquire images 
of the field of SN 2003gd at 24 /im on 2005 January 23 and 26 (days 678 and 681). In our 
measurement and analysis of the days 496/9, 670, and 678/81 observations, we used the 
same data as were available to S06. SN 2003gd was again observed at 24 /im within the 
^^jitzer supernova program of Sugerman et al. (PID: 30494) on 2006 September 1 (day 1264). 
Subsequent to the initial submission of this paper. Dr. Ben Sugerman kindly made this image 
available to us. We therefore also consider the implications of this observation. 



3. Results 

A point source at the SN position is clearly visible in the day 496/9 image in all four 
IRAC channels, with strong fading by day 670. This is illustrated in Figure [T](a,b), where 
we show the 8 /xm IRAG images from day 496/9 and day 670. Given the large decline in 
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flux, we deduce that most of tlie point-source flux detected on day 496/9 was due to tlie 
SN. However, measurement of tlie SN flux is challenging owing to the bright, complex field 
within which it lies. 



3.1. Day 496/9 Results 

3.1.1. PSF-Fitting Measurements 

We used the SNOOPY point-spread function (PSF) fitting package to determine the 
SN fluxes and coordinates. SNOOPY was originally designed by F. Patat to carry out SN 
photometry. It was implemented in IRAF by E. Cappellaro and is based on DAOPHOT, 
and has been tested and improved over a number of years. Several suitable PSF stars are 
selected in order to build the model PSF and measure the full-width at half-maximum in- 
tensity (FWHM). First a polynomial surface, of orders between 3 and 6 in x and y, is fitted 
to the background in a (10 x FWHM) x (10 x FWHM) region centered on the SN position, 
excluding the innermost square region around the SN, of side ~ 1.5 x FWHM. This is then 
subtracted from the image. Next the PSF fitting is performed on the SN. The fitted PSF 
is subsequently subtracted from the data to produce a residual image. This is inspected by 
eye, and the fitting procedure repeated until a residual image is obtained where there is little 
sign of the original point source. The code returns the x and y position and the flux within 
the PSF. It also provides a statistical uncertainty which is a measure of how well the model 
PSF describes the flux value and distribution at the SN position. However, the flux values 
are quite sensitive to the fitting of the image background with the polynomial surface. This 
may introduce an additional uncertainty in the absolute flux values, although the effect on 
the shape of the spectral energy distribution is likely to be less than this. 

The results and estimated uncertainties are shown in Table [TJ As a check of the PSF- 
fitting procedure, field stars were also measured using both this method and aperture pho- 
tometry. The aperture radius was 10" with a 15-20" concentric sky annulus. No significant 
systematic flux difference was found between the two methods at any wavelength. The 
root-mean-square (rms) scatter in the differences was 0.05-0.18 mag. The rms scatter at 
each wavelength was adopted as the uncertainty. As a further check, we performed aperture 
photometry for three stars in both post-basic calibrated data (PBCD) and SINGS-processed 
IRAC 8 fim frames, and found that the photometry agrees to within 5%. This test was 
also applied to the day 678/81 MIPS 24 /xm frames (see below) and similar consistency was 
obtained. 
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3.1.2. Image- Subtraction Measurements 



We also determined the difference in the fluxes between day 496/9 and day 670 via 
image subtraction. While this only gives the change in flux between the two epochs, it is 
a particularly powerful method of removing the effects of a spatiall y varying background 
such as is encountered in SN images from Spitzer flMeikle et al. II2006I ). Also, given the very 
weak flux at the SN location on day 670, this procedure provides a robust check on the net 
supernova emission. 



For each channel, the day 670 image (PBCD processed) was subtracted from the earlier 
Enhanced Data Product SINGS image through the use of image mat ch ing and subtr act ion 
techniques as implemented in the ISIS 2.2 image -subtraction p a ckage (lAlard 1 120001 ) . The 



8.0 fim subtracted image is shown in Fig. Wic). In lMeikle et al. I (120061 ) we demonstrate the 



applicability of the image-subtraction technique for Spitzer /IR AC SN data and investigate 
its uncertainties. Aperture ph otometry of the subtracted im ages was then carried out us- 
ing the Starlink package GAIA (jPraper. Gray, fc Berry II2002I ). A circular aperture of radius 
2'.'25 was used for the photometry. This aperture was chosen as a compromise between max- 
imizing the sampled fraction of source flux (the radius of the flrst diffraction minimum at 
the extreme red end of the 8.0 /xm channel is ~ 2'.'6) and minimizing any extended residual 
emission in the subtracted image. Aperture corrections were derived from the IRAC PSF 
images available on the Spitzer website. The correction factors were 1.23, 1.26, 1.50, and 
1.65 for 3.6, 4.5, 5.8, and 8.0 fim, respectively. 



For each measurement, the aperture was centered on the SN image using a combination 
of centroid estimates and visual inspection. The residual background level was measured 
using a clipped mean sky estimator and a concentric sky annulus having respective inner 
and outer radii of 1.5 and 2.8 times the aperture radius. The results are shown in Table [H 
The uncertainty was determined from the sky variance within the sky annulus. These error 
estimates were conflrmed by measuring the variance in the (day 496 - day 499) subtracted 
frame for each band, assuming a similar underlying error in the unsubtracted day 670 frame, 
and appropriately combining the two errors. These uncertainties are quoted in Table [TJ 
However, it is likely that additional systematic errors were present due to image-matching 
uncertainties. 



At 3.6 fim, 4.5 /xm, and 5.8 /xm the flux differences between the two methods all have 
less than 3a signiflcance (see Tabled]). At 8.0 /im the difference is over 4(T. As discussed 
below, we attribute this signiflcant difference to the presence of a residual source in the 



- 8 - 



day 670 image. In Tabled] we also show the PSF-derived IRAC fluxes obtained by S06 for 
day 496/9. There is reasonable consistency with our PSF results, although at 8.0 fim we 
see a higher flux at a significance of just under 3cr. Given the complexity of the field this 
difference is, perhaps, not too surprising. 

3.2. Days 670, 678/81 Results 

On day 670, there was no detectable source at or near the SN position at 3.6 fim, 4.5 fim, 
or 5.8 fim. However, sources were detected near the SN position at 8.0 fim on day 670 and 
24 fim on day 678/81 (see Figured]). We compared the positions of these sources with that of 
the SN. The coordinates of the SN were measured by applying PSF fitting (using SNOOPY) 
to the subtracted 8 /xm image (which was in the coordinate system of the day 496/9 SINGS 
image). For this, we used the PSF obtained from the day 496/9 8 /^m (SINGS) image which 
has a PSF identical to that of the subtracted image. We also measured the SN coordinates 
in the same image using three other methods: centroiding, optimal filtering, and Gaussian 
fitting as implemented in the IRAF CENTER task. The mean and standard deviation of 
the results from these four methods were adopted as the SN position and uncertainty, re- 
spectively. 

To convert the SN coordinates to the day 670 IRAC and day 678/81 M/PS" images, we 
derived geometric transformations between these images and the day 496/9 8 SINGS 
image. The transformation between the day 496/9 and day 670 8 fim images was obtained 
using the centroid coordinates of 20 isolated sources visible in both frames. We used IRAF 
GEOMAP to derive a general transformation including shifts, scales, and rotations in x and 
y, and a second-order polynomial for the nonlinear part. The transformation between the 
day 496/9 8 nm SINGS image and the 24 /xm MIPS image was obtained in a similar manner, 
using the centroid coordinates of 20 isolated sources visible in both frames. 

To measure the coordinates of the 8 fim and 24 /xm sources detected near the SN po- 
sition, we again used PSF fitting (SNOOPY). The associated uncertainties were estimated 
by simulating point sources using a Gaussian PSF with a flux similar to the faint source. 
Artificial sources were placed in each of the images at nine positions where the background 
was judged to have a similar level and complexity to that of the SN location. The coor- 
dinates of these sources were measured with PSF fitting using the same polynomial orders 
for modelling the background as for the actual 8 /xm and 24 fim sources. Finally, the mea- 
sured coordinates were compared with the known positions of the simulated sources and 
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the standard deviations of their offsets were adopted as tfie uncertainties in our coordinate 
measurements. 

Our conclusion from the above astrometric measurements is that the 8 /im source co- 
incides with the SN position to within 1" (90% confidence), and the 24 //m source coincides 
with the SN position to within 2" (90% confidence) i.e. in both cases the coincidence is to 
within 1 native pixel (1'.'2 at 8 /im, 2'.'5 at 24 /xm). The bulk of the position coincidence 
uncertainty arose from the PSF fitting to the days 670-81 sources with a smaller contribu- 
tion from the co-ordinate transformation and a negligible contribution from measuring the 
position of the SN on day 496/9. 

We measured the fluxes of the days 670-81 sources using our PSF-fltting procedure, 
obtaining 73 ± 7 /xJy at 8.0 /xm and 380 ± 90 /iJj at 24.0 /im. The 8.0 fim flux is consistent 
with the 77 ± 17 /xJy difference in the flux between the day 496/9 image and the subtracted 
image (see Table 1), indicating that the difference was due to the residual source in the 
day 670 image. The sensitivity on day 670 and day 678/81 is dominated by the effects of 
the bright nearby sources on the PSF fltting. For the other three IRAC channels, upper 
limits were obtained based on direct PSF measurements of the day 670 images and on the 
difference between the day 496/9 images and subtracted images. Our 2a upper limits at 
3.6 /im, 4.5 /xm, and 5.8 /xm are, respectively, 10 /^Jy, 20 /xJy, and 35 /xJy, rounded to the 
nearest 5 /xJy. 

3.3. Day 1264 result 

As indicated in Section 2, we were recently given access to the PBCD 24 yum MIPS 
image of SN 2003gd obtained within the Sugerman et al. S'piteer program (PID: 30494) on 
day 1264. Visual inspection suggests that the source near the SN position had faded since 
day 678/81. To investigate this more quantitatively, we subtracted the day 1264 image from 
the day 678/81 data using the procedures described in §3.1.2. 

We performed the subtraction on both the S'/A^GS'-processed day 678/81 image (pixel 
size 1'.'5) and on the two original PBCD images with the native pixel scale of 2'.'5 (an average 
of the native scale subtracted images was formed). A discrete source close to the SN position 
was observed in both subtractions. Aperture photometry of the source was carried out using 
a 6'/l radius aperture. The background was determined by using concentric sky annuli in the 
ratio 1.5:2 of the aperture radius, and also by placing the aperture (without sky annuli) at 



- 10 - 



a number of positions in a 2' x 0.5' box centered on the SN. The uncertainty was estimated 
from the rms value of the aperture values in the second method. The whole procedure was 
then repeated with a 4'.'4 aperture radius. Generally consistent results were obtained. The 
mean flux measured was 295 ± 70 /iJy. We conclude that it seems likely that the day 678/81 
24 fim source faded significantly by day 1264. 



4. Analysis 



The mid-IR fluxes for SN 2003gd at day 496/9 (Table OP) are plotted in Figure H The 
crossbars give the IRAC filter bands and the flux error bars are la. We show both the 
PSF-derived points a nd those derived by image subtraction. Also shown are BVRI points 
obtained on day 493 Jflendrv et al. Il2005h adjusted to day 49 6/9 using the SN 1987A hght 



curves . All the SN 2003gd points were dereddened using the ICardelli. Clavton. fc Mathis 



(|l989h extinction law with Ry = 3.1 and E{B - V) = 0.135 mag Jflendrv et al. Iboosj )^ 



There is clearly a strong mid-IR excess. The IR excess might be produced by an IR echo 
from circumstellar dust, but S06 argue that the decline rate is too high to be a typical IR 
echo. We find that it is, in f act, possible to reproduce the decline rate using a simple IR echo 
model jMeikle et al. I hood ) with a modest dust shell, although the shell parameters have 
to lie within quite a narrow range. Without more extensive temporal coverage, it is not 
possible to conclusively eliminate a significant IR echo contribution to the mid-IR emis- 



sion at day 496/9. Ho wever, the observed optical attenuation effects ( iHendry et al. 112005 



Sugerman et al. 1 120061 ) show that some dust condensation in the ejecta must have taken 



place. I n addition, there is no sign of ra dio emission, implying a paucity of circumstellar 
matter (IVan Dyk. Li. &: Filippenko II2003I ). Consideration of the deposited radioactive en- 
ergy also tends to support dust condensation at this epoch (see below). Given these facts, 
plus the need for a rather specific GSM shell geometry for an IR echo to reproduce the de- 
cline rate between day 496/9 and day 670, we proceed on the assumption that the day 496/9 
mid-IR flux was probably dominated by emission from newly formed dust in the ejecta. 



4.1. Comparison of SN 2003gd on Day 496/9 with SN 1987A 



In order to interpret further the day 496/9 mid-IR emission from SN 2003gd, ideally 
we would compare its spectral energy distribution (SED) with similar-epoch spectra from 
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a sample of SNe II-P, but such a database covering the 3-9 /xm range does not yet exist. 
The only pre-Spitzer 3-9 /im CCSN spectra are for SN 1987A. While SN 1987A was ini- 
tially atypical (i t arose from a blue supergiant star), its nebular optical/near-IR behavior 
has been shown (jPozzo et al. II2006I ) to be more similar to that of a normal SN II-P such as 
SN 2002hh. In addition, quite similar mid-IR spectral behavior ha s been found for SN 1987A 
and the Type II-P SNe 2004dj and 2005af around days 200-250 jKotak et al. Il2006h . Only 
in the [Ar III] 6.99 /xm line is significantly different behavior detected. We conclude that the 
nebular mid-IR behavior of SN 1987A is similar to that of Type II-P events like SN 2003gd. 



For comparison with SN 2003gd , we used SN 1987A spectra at 0.3-1.1 /im [SUS- 



( iBouchet fc Danziger 



1993 



PECT database andlPun et al. I fll995h]. 1.05-4.1 urn fMeikle et al. Ill993h. and 4 3-13.0 jjm 



Roche. Aitken. fc Smith 



1993; 



Wooden et al. Ill993l ). We used 



data from SN 1987A epochs as follows: optical/day 498, nea r-IR/day 494, mid-IR/day 494 



(iMeikle et al. Ill993l ). day 517 (IRoche. Aitken. fc Smith Ill993l ). and day 518 (IBouchet fc Danziger 



19931 ). This still left the blue half of the 8.0 fim band unrepresented. To fill in this gap. 



which includes the str ong [Ni II] 6.63 fim line, we used the day 415 SN 1987A KAO spectrum 
(IWooden et al. 1119931 ) . scaled and shifted to match the day 517/8 spectra in the overlap re- 
gions (4.5-5.3, 7.8-12 /xm). All the spectra were dereddened. In addition, to convert the 
spectra to the SN 2003gd epochs, small scaling adjustments were made using the SN 1987A 
light curves. The SN 1987A spectra were then scaled by 2.9 x 10~^ and 0.21 to compensate, 
respectively, for the distance and ^^Ni mass differences between the two SNe. 



Following all these adjustments, we found that the optical spectrum showed good con- 
sistency with the SN 2003gd photometry at both epochs. However, to match the SN 2003gd 
fluxes in the region of the IR excess (3-13 /im), we had to further increase the SN 1987A 
spectral fluxes by a factor of 2.0 for the image-subtracted points, and by a factor of 3.1 for 
the PSF-derived points. A compromise factor of 2.8 was applied (Fig. [2]). 



Comparison of the coeval SN 1987A IR spectrum with the day 496/9 SN 2003gd pho- 
tometry (Fig. [2D shows that much of the IRAC fluxes are likely to be due to emission from 
CO (~4.8 /im), SiO (~8.2 /im), fine-stru cture lines, and Br a . In the Kuiper Airborne 
Observatory (KAO) study of SN 1987A JWooden et anil993h . barely 20% of the fluxes 
corresponding to the IRAC 4.5 /xm and 8.0 /im bands were ascribed to emission from ejecta 
dust. On the other hand, in SN 2003gd the 4.5 /im point is not as far above the continuum 
as one might expect given the level of the SN 1987A CO emission. A similar, but less pro- 
nounced, effect may be apparent at 8 /im. This suggests that, while the factor of 2.8 scahng 
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is appropriate to match the IR continua of SN 1987A to that of SN 2003gd, it exaggerates 
the hne and molecular emission from the latter SN. Nevertheless, the non-dust contributions 
to the 4.5 fim and 8.0 fim fluxes of SN 2003gd are probably still significant, so these points 
should not be used for matching any dust emission model. In contrast, the 3.6 /xm and 
5.8 /im fluxes lie quite close to the SN 1987A continuum. Given that this continuum was 
due to emission from ejecta dust, we conclude that these points provide a fair measure of 
the emission from newly formed dust in SN 2003gd. We make use of these two points in 
matching the dust emission model. 



4.2. Dust Mass at Day 496/9 

To estimate the dust mass produced in SN 2003gd, we compared a simple analytical IR- 
emission model with the observed SEDs. An additional component was added to represent 
continuum emission from hot, optically thick gas in the ejecta. To select the likely grain den- 
sity distribution and grain materials for the dust emission model, we sought guidance from 
dust condensation calculations and the explosion models upon which they are based. Only 
a few papers have been published which describe SN dust condensation based on explosion 
models. Such papers fall into two categories: SN 1987A, and high-redshift low-metallicity 
progenitor SNe. No calculations for local Type II-P events have been published. We judge 
the SN 1987A dust models as probably being the more relevant. 



Kozasa. Hasegawa. fc Nomoto I (119891 ) and lTodini fc Ferraral (l200l[ ) have calculated dust 



condensation within the ejecta of SN 1987A. Th ese authors used the ejecta chemical com- 



position as determine d in nucleosynthesis models (IHashimoto. Nomoto. fc Shigeyama Ill989 



Nomoto et al. Ill99ll ). Both sets of authors adopted complete chemical mixing within the 



dust- 
tion. 



brming zone. Within this zone Todini & Ferrara assumed a uniform density distribu- 



Kozasa. Hasegawa. fc Nomoto (119891) used the density profile from an explosion model 



( IHashimoto. Nomoto. fc Shigeyama I Il989l ) but this also is roughly flat. Similar dust-type 
abundances were obtained by both sets of authors, but neither make explicit predictions 
about t he dust distribution w ithin the ejecta. Recent three-dimensional CCSN explosion 
models (IKifonidis et al. II2006I ) confirm that extensive mixing of the core takes place. The 
same models also show that the density structure is likely to be exceedingly complex, with 
high-density clumps moving out through lower-density gas. How this affects the dust distri- 
bution has yet to be determined. 



Given the current state of knowledge, we assume that dust of uniform number den- 
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sity forms throughout the zone containing abundant refractory elements. The extent of 
this zone can be assessed us ing the late-time wi dths of metal lines. In the day 493 opti- 
cal spectrum of SN 2003gd (iHendry et al. II2005I ). the maximum velocities implied by the 
metal lines generally do not exceed ~2000 km s~^. This upper limit is adopted as the 
size of the dust-forming region. The uniform density assumption is conservative in that 
it provides the least effective way o f hiding dust grains in optically thick regions. Guide d 



by the dust-formation calc ulations (IKozasa. Hasegawa. fc Nomoto I Il989l : iTodini fc Ferrara 



2OOII : iNozawa et al. II2003I ). we included silicate, amorphous carbon, and magnetite dust 
in the mass ratios 0.68/0.16/0 .16. The mass absorption functions for the three materials 



were taken from the literature (iLaor fc Draine Ill993l : iRouleau &: Martin Ill99ll : iKoike et al. 
198lh . 



Our dust IR-emission model compr ises a uniform sphere of isothermal d ust grains. Fol- 
lowing the escape probability formalism ( Osterbrock Ill989 : Lucy et al. Il989h . the luminosity 
{L^) of the sphere at frequency z/ is given by 



U = 2n'R'B,{T)[T-\2T'^ - 1 + (2r, + 1)6-^-)], 



(1) 



where R is the radius of the dust sphere at some time after the explosion, By{T) is the Planck 
function at temperature T, and t^, is the optical depth to the center at frequency v. For a 
grain size distribution dn = ka~"^da, where dn is the number density of grains having radius 
a a + da, m is typically between 2 and 4, and k is the grain number density scaling factor, 
it can be shown that Ti, = ^TikpKyRj^[a^^ — a^^], where p and are, respectively, the 
density and mass absorpt ion coefficient of the grain material . The grain size distribution 
law was set at m = 3.5 (IMathis. Rumpl. fc Nordsieck 1 119771 ) with = 0.005 /im and 
Q mrt.r = 0.05 Atm. The total mass of dust, M^, was then found from Md = AtiR'^t^/Sku 



flLucv et al. \\l98% 



The model-free parameters are the grain temperature, sphere radius, and grain number 
density scaling factor, k . Th ese were adjusted to reproduce just the 3.6 pm and 5.8 pm 
points. IWooden et al. I (119931 ) showed that during the second year of SN 1987A, the dust- 
emission continuum could be contaminated by blackbody emission from hot, optically thick 
gas, as well as by free-bound radiation. Here we represent both effects using a single hot 
blackbody. The hot component was adjusted to match the underlying continuum of the 
scaled SN 1987A optical spectrum and not the broad-band points of SN 2003gd, which 
would contain a significant contribution from the many strong emission lines. We found that 
the effect of the hot component on emission longward of 3 pm was small. 
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Model matches to both the PSF-fitting and image-subtraction-derived fluxes were ob- 
tained. The dust emission models are shown in Figure [2l We found that to achieve reasonable 
matches to the data it was necessary to increase the dust mass until it was optically thick in 
the mid-IR. Consequently, we were unable to derive a unique dust mass since, as we increase 
the optical depth, ever larger amounts of dust can be "hidden" with little effect on the ob- 
served radiation. We therefore, conservatively, sought the minimum dust mass which would 
provide a satisfactory match to the data. The model parameters including the derived dust 
masses are given in Table [H 



Dust masses of 6 x 10~^ Mq (PSF fitting, day 496/9 image) and 4 x 10"^ Mq (aper- 
ture photometry, subtracted image) were obtained. A 1 Mpc reduction in distance reduces 
the masses by about 10%. The uniform dust distribution of our model, optically thick at 
10 /im, would surely extinguish all r netal lines in the optica l region. Yet, as late as day 493 
(IHendry et al. I l2005l ) and day 521 (ISugerman et al. II2006I ). such lines could still be seen. 
This implies that the dust distribution must have been "clumpy," allowing some of the op- 
tical line radiation to escape from the nebula. The presence of clumping is confirmed by 
consideration of the "covering factor," /. This is obtained by dividing the projected area 
of the model dust sphere by the projected area corresponding to the estimated extent of 
the dust-forming zone (2000 km s~^). A covering factor of ~0.15 was ob tained (Table El). 



This may also account for the relatively modest extinctions in the R band (ISugerman et al. 



2006 ). We note that SN 1987A showed strong evidence for dust clumping (ILucy et al. 



199ll ). In Section 5, we suggest that, in general, SN ejecta dust becomes optically thick in 
the mid-IR when the dust mass exceeds only a few times 10"'^ Mq. 



It is argued above that the 8 /im point should not be used for constraining the dust 
model due to possible contamination by SiO emission. Nevertheless, we investigated the 
effect of including this point and found that satisfactory matches can be obtained using 
somewhat lower temperatures and higher radii for the model. Similar dust masses are de- 
rived. However, the match to the SN 1987A spectrum redward of 8 /im is very poor, with 
the model flux exceeding the continuum flux by about a factor of two. Given the argument 
in § 4.1 that the nebular behavior of SNe 1987A and 2003gd is similar, we conclude that SiO 
is indeed contaminating the SN 2003gd spectrum. 

The total luminosity of our dust model (for the match to the PSF-derived fluxes) 
plus the estimated total optical/near-IR contribution (i.e. line/molecular emission plus 
underlying continuum) is 2.1 x 10^^ erg s^^, with roughly 30% of the luminosity arising 
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from the dust. The ^^Ni mass inferred bv iHendrv et al. I (120051 ) is 0.016lo;oo6 M©. Divid- 
i ng the observed total luminos ity by the radioactive decay energy deposited in the ejecta 
fiLi. McCrav. fc Sunvaev I llQQsh . scaled to the ^^Ni mass of SN 2003gd, we obtain 1.2;°;^. 
Thus, the total luminosity is similar to that resulting from the deposited radioactive energy. 
This tends to support the proposition that newly condensed ejecta dust was responsible for 
the mid-IR continuum emission. Use of a lower distance would reduce the radioactive decay 
energy required to produce the observed flux. However, this would not significantly affect 
the energy constraints as the inferred ^^Ni mass would also fall — that is, the fraction of 
radioactive decay luminosity required to produce the observed flux would stay about the 
same. The dust masses we derive are about 25% of the 2.0 x 10~^ M© which S06 obtain 
from their smooth model fit at the same epoch (day 496/9). The S06 models are shown in 
Fig. [2] (dashed lines). Between 3 /im and 10 fim their models are in approximate agreement 
with ours, but at longer wavelengths our model shows a much sharper decline. It appears 
that the S06 model predicts a component of colder dust and this would account for their 
larger dust masses. Their model invokes a source luminosity of 2.6 x 10^^ erg s~^. This is 
~20% larger than the value in our model but is still consistent with the radioactive energy 
deposited, given the uncertainties in the ^^Ni mass. 



4.3. The Days 670-681 Sources 

S06 found a 24 /im flux on day 678/81 of 106 ± 16 /xJy. They also reported upper limits 
at 3.6 /im, 4.5 /im, 5.8 /im, and 8.0 /xm on day 670. It is from these later-epoch measure- 
ments that they deduce a dust mass of 0.02 Mq. As indicated above, we also obtained no 
detection at 3.6 /xm, 4.5 /xm, and 5.8 /im on day 670. However, at 8.0 /im we obtained a 
significant detection in our PSF fitting of 73 ± 7 /xJy. Moreover, inspection of Fig. [Hb) does 
appear to confirm the presence of a source close to the SN position. At 24 /im our mea- 
sured flux of 380 ±90 /iJy (see above) is about a factor of 4 larger than that obtained by S06. 

To investigate this flux difference we assessed the day 678/81 MIPS sensitivity at the 
source position using a number of methods. The complex field in the SN vicinity makes direct 
noise estimation quite difficult. Therefore, to determine the underlying pixel-to-pixel noise, 
we subtracted the day 681 image from the day 678 image using the procedures described 
above. We then measured the noise at the SN location. We used a 6'.'1 radius aperture which 
encompasses about 0.93 of the flux in the Airy disk at 24 /im. The flux in the aperture 
was measured at a series of locations within 40" of the SN position. The rms value, after 
aperture correction, is ~200 /iJy. However, as the subtracted image contained the noise of 
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the two original images, we divided this by y^, yielding 140 /iJy as the intrinsic sensitivity 
(la) of the MIPS data. A further a/2 improvement of the sensitivity to 100 fiJy arises from 
the fact that the SINGS-processed image, used by S06 and ourselves for the PSF fitting, is 
a combination of the two SINGS PBCD images. 

As a further check, we examined the MIPS sensitivity in the SINGS-processed image 
well away from the galaxy in a relatively "clean" part of the sky, lying about 4' south of the 
SN location. Artificial stars were placed at 11 different positions within a 1' x 3' area. The 
input star flux was set at 130 /iJy and the flux at each of the 11 positions was measured 
by aperture photometry, using a 6'.'1 radius aperture and a sky annulus between 1.5 and 
2 times this radius. The effective sensitivity was assessed from the dispersion in the flux 
values. The procedure was repeated with an input star flux of 5300 /iJy. The dispersion 
in both the low-fiux and high- flux cases gave about the same result, indicating that even 
well away from the galaxy, the sensitivity was background limited. The la sensitivity, after 
aperture correction, was found to be about 90 /iJy, similar to the value obtained from image 
subtraction. 



As a final check on the above procedures, we used the Spitzer Sensitivity-Performance 
Estimation Tool (PET) to estimate the intrinsic sensitivity. The measured background near 
the SN is about equivalent to the "high background" setting of the PET. From this we 
derive an intrinsic la sensitivity of 60 /iJy, of similar magnitude to the directly-determined 
sensitivity values. 

We conclude that the actual la sensitivity of the day 678/81 MIPS image at the SN po- 
sition was ~90 /iJy, much larger than the 16 /iJy claimed by S06. Moreover, the 106 /iJy flux 
at the SN position claimed by S06 would only yield a S/N of about unity, not ~ 6.5 as they 
reported. However, we note that scaling the S06 result by a factor of 4 yields 424 ± 64 /iJy. 
This is more consistent with both our flux value and with our separately measured MIPS 
sensitivity. We suspect, therefore, that there is an error in the 24 fim flux reported by S06. 



What is the origin of the 8 /xm and 24 /im sources on days 670-81? Given the very 
crowded field within which the supernova occurred, and the fact that CCSNe tend to occur 
near star-forming regions, a cool background source lying close to the SN might be con- 
sidered. A similar situation wa s described recently for Spitzer observations of the CCSN 
SN 2002hh jMeikle et al. Il2006h . However, the fading of the 24 source (§ 3.3) tends to 
rule out a background source, at least for most of the 24 fim flux. We shall therefore consider 



the implications of assuming that the sources are ultimately due to the SN. 



We first hypothesize that the days 670-81 8 fim and (unsubtracted) 24 fim sources have 
the same origin and that this origin is the SN ejecta. A simple blackbody match to our 
days 670-81 flux measurements yields a temperature of 250 K, a radius of 1.9 x 10^^ cm, and 
a luminosity of 9.7 x 10^^ erg s~^. This is immediately problematic. To attain a radius of 
1.9 X 10^^ cm the material at the outer limit of the blackbody would have to be travelling 
at 3200 km s~^. This is substantially larger than the 2000 km s~^ limit on metal velocities 
indicated by late-time spectra. In addition, after adding an additional 1.1 x 10^^ erg s~^ due 
to the optical/near-IR emission estimated from the optical photometry, we obtain a total 
luminosity of 10.8 x 10^^ erg s~^. This is a factor of 4 more than the total radioactiv e deca y 
energy deposited in the ejecta, according to the formula of iLi. McCray. fc Sunyaev I (119931 ). 
scaled to the ^^Ni mass of SN 2003gd. Indeed, it exceeds the tota/ radioactive luminosity (i.e., 
including escaping gamma rays) by more than a factor of 2. Even allowing for the uncertainty 
in the ^^Ni mass, the energy deficit is severe. (We note that, even with their apparently un- 
derestimated 24 fim flux, the day 678/81 luminosity invoked by S06 exceeds the deposited 
radioactive energy by ~50%.) It is possible for the bolometric luminosity to exceed that of 
the instantaneous radioactive decay deposition when the recombination timescale exceeds 
the radioactive or expansion timescales. However, this cornmenc es at much later epochs 
(>day 800) than those considered here (IKozma fc Fransson Ill998l ). Thus, on both energy 
and velocity considerations, we have evidence that most of the 8 fim and 24 /xm fluxes cannot 
be due to emission from supernova ejecta dust. 



Let us now suppose that only the fading component of the 24 /im source is due to 
ejecta dust while the remainder of the 24 fim flux plus some or all of the 8 fim flux is 
due to a background source. At 250 K, to match the fading component would require a 
blackbody luminosity of 3 times the likely deposited radioactive luminosity and a veloc- 
ity of 2800^^1^1^ km s-\ where the error is due to the flux uncertainty. Even if we reduce 
the distance by 1 Mpc and use the lower limit of the flux values, the velocity still exceeds 
2000 km s~^ and the luminosity still exceeds the deposited radioactive luminosity by a factor 
of 2. Reducing the temperature from 250 K to 150 K, the luminosity falls by 30% but the 
velocity of the blackbody surface rises to an increasing ly implausible 6700l^|^|] km s-\ (We 
note that, in their model, S06 invoke an outer limit for their dust zone of ~8000 km s~^, 
which is even more unlikely.) Increasing the temperature above 250 K also does not help 
since the luminosity deficit problem would worsen. Moreover, this would produce an 8 fim 
flux in excess of that seen near the SN position. Similar results are obtained if we employ 
our dust emission model rather than a blackbody. We conclude that most of the fading 
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component of the 24 /im flux cannot be due to dust in the SN ejecta. 



If the mid-IR flux observed near the position of SN 2003gd is not due to condensing 
dust in the ejecta, then what could be the origin of the emission? The substantial fading 
at 24 fim between day 678/81 and day 1264 points to a causal connection with the SN. A 
possible scenario is that the mid-IR emission originated in an IR echo from circumstellar or 
interstellar gas. As an illustration, we have estimated the parameters of a dust sheet lying 
in front of the SN required to reproduce the fa ding component of t he 24 /xm flux. We used 
an IR echo model similar to that described by Meikle et al. (I2OO6). T he input bolometric 
light curve was based on the information given by Hendry et al. I ( 2005 ). with a single grain 
radius of 0.07 fim. Estimates were repeated using the specific grain emissivities of different 
dust species. Preliminary results suggest that the 24 /im flux can be reproduced with a dust 
sheet of H number density ~10 cm~^, gas-to-dust ratio of 100, lying 10-20 pc in front of 
the SN. At this distance the dust temperature is 75-90 K. The optical depth to UV-optical 
photons is ~0.2. The echo radius would be about O'.'l and so such a source would be effec- 
tively coincident with the SN position. To account for the fading the dust sheet would have 
to be of irr egular density on scales of a few p arsecs (a fraction of an arcsecond). ISugerman 
( I2OO5I ) and IVan Dyk. Li. fc Filippenko I (120061 ) found an optical echo on day 623 lying at O'.'S 
from SN 2003gd, with a strong concentration to the NW. They showed that this could be 
explained by a dust sheet lying about 100 pc in front of the SN. The one-sided appearance of 
the optical echo suggests that such dust sheets can indeed exhibit large density fluctuations 
on a scale of only a few tenths of an arcsecond. We conclude that an IR echo may well be 
responsible for the variable component of the 24 /xm flux from SN 2003gd. Further discussion 
of the IR echo hypothesis as applied to SN 2003gd is beyond the scope of this paper. 



The key point following from the above discussion is that most of the mid-IR flux at 
days 670-681 could not have been due to dust in the supernova ejecta. In particular, it 
suggests that the inference by S06 of a large mass (0.02 Mq) of ejecta dust is unjustified. 



5. Conclusions 



We have examined late-time mid-IR observations of the Type II-P SN 2003gd and find 
the following. 
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(i) By day 496/9, at least 4 x 10"^ M© of dust had formed in the ejecta of SN 2003gd. 
The larger (factor of ~4) mass indicated by the smooth model of S06 appears to be due to 
the presence of a larger component of cold dust, but this has no direct observational support. 
After allowing for differences in ^^Ni production, we find that the optical flux of SN 2003gd is 
similar to the coeval value for SN 1987A, while the 3-9 /xm flux is almost three times stronger. 
This may indicate more efficient dust production in SN 2003gd. Nevertheless, there is no ev- 
idence at this epoch that the absolute dust production was u nusually high. The du st masses 



and temperatures are similar to those inferred for SN 1987A (jWooden et al. Ill993l ). There is 
also evidence that the dust in SN 2003gd formed in clumps. Comparison with coeval spectra 
of SN 1987A shows that even as late as day ~500 the extraction of information about dust 
formation from broad-band photometry has to be approached with caution due to the effects 
of other emission mechanisms. This underlines the desirability of acquiring low-resolution 
spectra for such studies, since this would allow correction for forbidden lines and molecular 
emission. 



(a) Emission from point-like sources close to the SN position was detected on day 670 
at 8 /im and day 678/81 at 24 /xm. The fading of the 24 /im source (§ 3.3) tends to rule out 
a background source as the origin of the mid-IR emission, at least for most of the 24 /xm 
flux. However, energy and velocity considerations also rule out the ejecta of SN 2003gd as 
the origin of most of the mid-IR fluxes. The inference by S06 of 0.02 of ejecta dust is 
based on a 24 /xm flux which we find is only a quarter to a third of the true value. But even 
if we adopt their flux, their claim of such a large mass of dust is unconvincing. The large 
dust mass they flnd appears to be a consequence of the low characteristic temperature in 
their model. However, in order that sufficient mid-IR radiation should escape, it seems that 
the dust formation zone has to be as large as ~8000 km s~^, in conflict with the observed 
metal line velocities. Also, in spite of the low temperature, the input luminosity of their 
model required to reproduce the 24 /xm flux exceeds that of the likely deposited radioactive 
decay luminosity. Had the correct (much larger) 24 /im flux been used, these difficulties 
would have been even greater. We conclude that the 0.02 of ejecta dust deduced by S06 
is unsupported by the data. These Spitzer observations provide no basis for the S06 claim 
that "the [dust] condensation efficiency implied by SN 2003gd is close to the value of 0.2 
needed for SNe to account for the dust content of high-redshift galaxies." There is, as yet, 
no direct evidence that CCSNe are major dust factories. 



An additional argument against a large detected mass of ejecta dust in SN 2003gd is 
pointed out by the referee. The mass of dust inferred by S06 for SN 2003gd is a factor of 25 
greater than the maximum amount in SN 1987A determined by lErcolano. Barlow, fc Sugerman 
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( 120071 ) using a similar model. These authors suggest that this implies a much higher conden- 
sation efficiency in SN 2003gd. Yet, the diminution of the optical light curve of SN 2003gd 
shown in S06 about 100 days after dust formation is only about 1.3 ma g, very similar to 



that for SN 1987A at about the same epoch after dust formation (e.g., IWhitelock et al 



19891 ). Also the blueward shifts of the emission lines are not greater in SN 2003gd than in 



SN 1987A (e.g.. iDanziger et al. Ill99ll ). Given the apparently very different estimated dust 
masses between the two SNe, it is difficult to see how these observed similarities are possible. 
One might argue that in SN 2003gd the dust is more concentrated to the center, but this is 
apparently belied by the shift in the Ha lines (S06) which presumably arise farther out in 
the envelope. 



The goal of determining the true dust production in SNe via the thermal emission from 
the grains is very challenging. Even at wavelengths as long as 24 /im, it is likely that the dust 
forming in SN ejecta would become optically thick long before a universally significant mass 
of dust was formed. For example, consider a u niform distribution o f astronomical silicate 
grains. For A > 20 /xm, Ki, ~ 1000(A(/im)/20)~^ (jLaor &: Draine Ill993l ). where n^y is the mass 
absorption coefficient (cm^ g~^) at frequency z/. If we set > 3, where is the optical 
depth to the dust sphere center at frequency u, and let the radius of the refractory element 
zone be as large as v = 3000 km s~^, we obtain from our uniform dust model 



Md > 1.5 X 10-2(A(/im)/20)2(t(days)/600)^ M©. 



(2) 



At 24 nm, and as late as 2 years after the explosion, the lower limit for the dust mass would 
still only be 3 x 10^'^ M0. Similar lower limits are obtained for other grain materials such as 
amorphous carbon. In general, rises toward shorter wavelengths, producing even smaller 
lower limits. Dust measurement at still later epochs becomes increasingly difficult as the 
grains cool beyond the sensitivity limits of even Spitzer. 



Of course, if the grains are arranged in optically thick clumps, then a large mass of 



in SN 1987A ( 


Lucv et al. 


1989 


Ercolano, Barlow, & Sugerman 



in just the optical/near-IR region, but thin in the mid-IR, then the dust mass may be con- 
strained by the observed luminosities in the two wavelength regions. However, once the dust 
becomes optically thick in the mid-IR it is possible to derive lower limits only. As explained 
above, this situation sets in when the dust mass exceeds only a few times 10"'^ M©, well 
below the cosmologically interesting limit of ~ 0.1 Mq. 
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The value of mid-IR studies of CCSNe, such as the Spitzer work described here, is that 
they can test whether dust condensation is common in typical events (i.e., SNe II-P). This 
is an essential step if we are to demonstrate that CCSNe are major sources of universal 
dust. But if large masses of dust are formed in SN ejecta, the direct measurement of the 
total masses involved is likely to require observations at much longer wavelengths. However, 
even with such observations there would remain the challenge of eliminating the effects of 
IR echoes. 
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A. SN Dust Production Required to Account for Observed High-Redshift 

Dust 



Models of dust formation in CCSNe (iTodini &: Ferrarall200ll : iNozawa et al. II2003I ) suc- 
ceed in producing copious amounts of dust — around 0.1-1 Mq even in the low-met allicity 
environments of the early universe. This corresponds to a SN dust condensation efficiency 
of about 0.2 (IMorgan fc Edmunds II2003I ). where the efficiency is defined as the dust mass 
divided by the total mass of refractory elements. This is enough to account for the quantity 
of dust seen at high redshifts. 



As a demonstration, let us consider the results of iBertoldi et al. I (120031 ). In their study 
of high-redshift quasars they deduced, from the far-IR luminosities, a star-formation rate of 
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~3000 M0 yr ^ and a dust-formation rate of ~1 yr ^. Consider a simple stellar mass 
spectrum 

dN = -fM-^-^dM, (Al) 

where dN is the number of stars in the mass interval M to M + dM, 7 is a constant with 
units of (mass)^'^, and all stars lie within the mass range 0.2 Mq < M < 30 Mq. The 
total stellar mass in the interval M to M + ciM is given by dMtot = MdN = 'jM-'^-^dM. 
Integrating this equation over the stellar mass range we obtain Mtot = 4.I7. In one year we 
have Mtot = 4.I7 = 3000 M0, so 7 = 732. Stars of mass exceeding about 8 M0 will end 
their lives as CCSNe. Integrating equation (Al) over the range 8 Mq to 30 Mq, we obtain 
= O.O38I7 = 28.0 — that is, about 28 CCSNe per year would occur. To produce 1 M© 
of dust per year, the average dust yield of eac h SN would have to be 0.036 M(7^. Thus, the 



production rates of dust condensation models (ITodini fc Ferrarall200ll : iNozawa et al. II2003I ) 



are indeed sufficient to account for the high-redshift dust. 



The progenitor of SN 2003gd had a mass in the range 6-12 Mq ( iVan Dyk. Li. &: Filippenko 



2003 



Hendry et al. This would produce about 0.3 Mq of refractory elements (IWoosley fc Weaver 



I995I ). Consequently, for such a SN to match the required average dust production, the re- 



fractory elements would have to be converted into dust with an efficiency of about 0.1. 
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Table 1: Mid-IR Photometry on Day 496/9 at the Position of SN 2003gd 



Author/method 




Flux 


(/^Jy) 






3.6 fim 


4.5 fim 


5.8 fim 


8.0 /xm 


This work: PSF fit. 


19.9(3.6) 


74(12) 


85.5(6.5) 


180(15) 


This work: Im. sub. 


15.6(1.3) 


72.2(2.3) 


60.2(7.1) 


103.2(7.7) 


S0(3: PSF fit. 


20.8(2.(3) 


73.8(5.6) 


61.9(7.3) 


103(22) 



Col. 1 gives the photometry method used to derive the fluxes. The "Im. Sub." values were obtained by 
aperture photometry of subtracted images (see text). Uncertainties {la) are shown in brackets. The fluxes 
obtained by S06 for this epoch are shown for comparison. 
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Table 2: Model Parameters for Day 496/9. 



Method 


^dust 




Dust Mass 


/ 


^hot 




(K) 




Mo 




(K) 


PSF fit. 


525 


2.6 


6 X 10-5 


0.17 


6700 


Im sub. 


525 


2.3 


4 X 10-5 


0.13 


6700 



Col. 1 gives the photometry method used to derive the fliixes to which the model was matched. Optical 
depths to the center at 10 ^m are shown in col. 3. The dust masses in col. 4 were derived assuming a distance 
of 9.3 Mpc. A 1 Mpc reduction in distance reduces the masses by about 10%. In col. 5, a "covering factor" 
/ is shown (see text). Col. 6 gives the adopted temperature of the hot component. 
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Fig. 1.— Panels (a) and (b) show the field of SN 2003gd at 8 fim (IRAC) on day 496/9 and 
day 670, respectively. Panel (c) shows the result of subtracting (b) from (a). The bright 
point source is SN 2003gd. Panel (d) shows the field of SN 2003gd at 24 /xm (MIPS) on 
day 678/81. North is up, and east is to the left. The pixel scales are Ql!75 pixel"^ and l'/5 
pixel"^ at 8 /im and 24 /im, respectively. Each panel is centred on the supernova position, 
determined using PSF fitting and the IRAF GEOMAP package (see text). 
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Fig. 2. — Spectral energy distribution of SN 2003gd on day 496/9 compared with the 
coeval SN 1987A spectrum and dust emission models. All data have been dereddened. The 
spectrum has been corrected for redshift. See text for explanation of interpolation of the 
data to days 496/9. Open (blue) squares in IR region: PSF-derived IRAC points. Solid 
(blue) squares: image-subtraction-derived IRAC points. The cross b ars indicate the IRA C 
bandwidths. Solid (black) squares in optical region: photometry from lHendry et al. I (120051). 
Faint structured (green) spectrum: optical/IR spectrum of SN 1987A scaled to allow for 
differences in distance and ^^Ni mass. The 3-13 /^m spectrum flux has been multiplied by 
an additional factor of 2.8. Light continuous line (red): dust emission model matched to 
PSF-derived fluxes at 3.6 and 5.8 fim (this work). Heavy continuous hne (red): match to 
image-subtraction-derived fluxes at 3.6 and 5.8 /im (this work). Note that the hot component 
was adjusted to match the underlying continuum of the scaled SN 1987A optical spectrum 
and not the broad-band points of SN 2003gd, which would contain a significant contribution 
from the many strong emission lines. Dashed lines (black): dust emission models from S06. 
Long-dash: smooth model. Short-dash: clumped model. 



